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Abstract. Directed polymers (strings) and semiflexible polymers (filaments) are 
one-dimensional objects governed by tension and bending energy, respectively. 
They undergo unbinding transitions in the presence of a short-range attractive 
potential. Using transfer matrix methods we establish a duality mapping for 
filaments and strings between the restricted partition sums in the absence and 
the presence of a short-range attraction. This allows us to obtain exact results 
for the critical exponents related to the unbinding transition, the transition point 
and transition order. 
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1. Introduction 

Directed polymers (or "strings" in the following) are one-dimensional objects governed 
by their tension which tends to minimize the contour length of the polymer. 
Semiflexible polymers (or "filaments" in the following), on the other hand, are 
governed by their bending energy which tends to straighten the polymer. In the 
presence of a short-range attractive potential, these objects undergo unbinding or 
desorption transitions which represent a number of important critical phenomena 
[TJH]- The unbinding of strings describes wetting 1 , polymer adsorption [3J, pinning 
of flux- lines in type-II superconductors [4], or roughening of crystal surfaces [5]. The 
unbinding of filaments describes adsorption and bundling of many biopolymers (DNA, 
F-actin, microtubules) and polyelectrolytes with large persistence lengths [6]. 

In this letter we use transfer matrix (TM) methods to derive a duality mapping 
for filaments and strings between the restricted partition sums in the absence and 
in the presence of a short-range attractive potential. This allows us to obtain the 
unbinding and desorption transition point, the order of the transition, and a set of 
scaling relations for the critical exponents of bound and unbound filaments and strings. 

2. Model 



We consider strings or filaments in 1 + d± dimensions which are oriented along the 
x-axis such that we can parameterize the contour by a dj_-dimensional field z(x) of 
displacements perpendicular to the x-axis with < x < L where L is the projected 
length of the string or filament. The Hamiltonian for strings is given by the sum 
of the tension energy f Q dx(a/2)(d x z) 2 with a string tension a and the potential 
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energy J Q V(z(x)), where V(z) contains an attractive potential well of range £ a which 
favours the configuration z = 0. The Hamiltonian for filaments is given by the sum 
of the bending energy dx{K/ 2){d x z) 2 and the potential energy V (z(x) , d x z) . 
k is the bending rigidity of the filament and L p — 2k /T the persistence length at 
temperature T. The expression for the bending energy in the parameterization by 
the projected length is appropriate if either the total length L or the longitudinal 
correlation length £|| to be defined below are small compared to L p . In contrast to 
the string, the filament has a well-defined tangent vector at each point, and therefore, 
also the external potential V(z, v) can depend on the tangent vector v = d x z. 

Generic potentials are of the form V = V r + V a + V p and contain a hard-core 
potential V r , a short-range attractive potential V a , and eventually a long-range power- 
law potential V p . The hard core potential V r is given by V r (z) — oo for |z| < £ r and 
V r (z) — otherwise. The short-range attractive potential V a has finite range £ a 
and a potential strength W < 0, i.e., V a (z) = W$(v) for |z| < £ a (£ a > £ r ) and 
V a (z) = otherwise. For strings we can only consider position-dependent potentials 
and set $(v) = 1. For filaments we include the dimensionless function $(v) modeling 
an additional orientation-dependence of the attractive potential. The potential V a 
attains the asymptotic form V a (z) = G£~ d± $(v)5(z) in the limit of small £ a where 
G = Wir d± / 2 /r(l + d±/2) < 0. Finally, we can also include attractive long-range 
power-law potentials V p (z) = w\z\~ p for |z| > £ a . Our results apply to potentials V p 
that decay sufficiently fast, i.e., potentials V p with p > 2 for strings and p > 2/3 for 
filaments [7]- 



3. Transfer matrix equations 

In order to simplify the notation, we introduce rescaled quantities measuring energies 
in units of the temperature T and lengths in units of T /2a for strings and in units 
of the persistence length L p = 2k JT for filaments. In rescaled units the restricted 
partition sum for strings with fixed initial point zo = z(0) and end point z = z(L) 
takes the form 

Zl{z\zq) = / T>z(x) exp < — / dx 
J(z -o) [ Jo 

In analogy with quantum mechanics, this path-integral fulfils a Schrodinger equation 
in imaginary time, the partial differential TM equation given by 

d L Z L = V\Z L - V(z)Z L (2) 

with the boundary condition Zq(z\zq) = 8{z — zq) at L = 0. The Laplace transform 
of the restricted partition sum with respect to L, Z s = J g dLe~ sL Zl, fulfils the 
differential TM equation 

sZ s - V 2 Z Z S - V{z)Z s + S(z - z ) (3) 

where the last term on the right hand side represents the boundary condition at L = 0. 
For a sufficiently attractive potential, there exist bound states for which we make the 
Ansatz Zl(z|zo) ~ i^e(^) exp(-EL) where E < is the free energy difference between 
the bound state and the free state (obtained for V = 0). The eigenfunction ^s(z) for 
the energy level E then solves the stationary Schrodinger equation 

- E^ E = Vl^ E - V(z)^e (4) 



-(d x zf + V(z(x)) 



(1) 
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with E < for a bound state. We impose the normalization J ip E ( z ) = 1. Then 
the solution satisfying the proper boundary condition is obtained by summing over 
all energy levels E n , Zl{z\z ) = J2n ^-E„ (z)^^ ( z o)e _B " i , where the ground state 
Eo dominates the sum for lengths L exceeding the correlation length £|| = l/l-Eol 
(assuming that binding is weak such that the continuous scattering spectrum starts 
at Ei = 0). 

For filaments we can proceed similarly starting from the restricted partition sum 
in rescaled units, in which we additionally fix initial tangent vo = d x z(0) and end 
tangent v = d x z(L). This partition function is given by 

p(z,v;L) ( pL 

Z L (z, v|z , v ) = / Vz(x) exp < - / dx 

J(z ,v ;Q) { JO 

and again fulfils a Schrodinger-like differential TM equation [8j [9] 

d L Z L = -v • V Z Z L + V 2 V Z L - V(z, w)Z L (6) 

with the boundary condition Zo(z, v|zo, vo) = S(z — zq)S(v — vo) at L = 0. As 
for strings we can consider the Laplace transform which fulfils the differential TM 
equation 

sZ s = -v • V Z Z S + V 2 V Z S - V(z, v)Z s + 5(z - z Q )5(v - v ) (7) 

where the last term on the right hand side stems from the boundary condition at 
L = 0. For sufficiently strong attractive potential, there exist bound states for which 
we make the Ansatz Zl(z, v|z , v ) ~ ^B(z,v)exp(-£'L), where E < is the free 
energy difference between bound and free state. The eigenfunction iPe(z,v) for the 
energy level E then solves the stationary version of the Schrodinger-like equation © 

- E^ E = -v • V z ip B + VlipE - V{z, v)^ E (8) 

with E < for a bound state. As for strings, we impose a normalization 
Jz Jv V'bC 2 ) v )V ; s(z, — v) = 1, and the solution satisfying the proper boundary 
condition is obtained by summing over all energy levels E n . For lengths L exceeding 
the correlation length £|| = I/I-EqI, the ground state dominates and Zl(z, v|zo, Vo) ~ 
ipE (z, v)V»b (z , -v )e- E ° L . 



-{dlzf + V{z{x)A^) 



(5) 



4. Scaling behaviour and exponents 

Strings and filaments differ in the scaling of free mean-square displacements, i.e., 
(|z| 2 ) ~ L 2| > for V — where £ is the roughness exponent. Strings show diffusive 
behaviour with £ = 1/2, whereas filaments have ( ~ 3/2. Tangent vector fluctuations 
scale as (|v| 2 ) ~ L 2 ^^ 1 ^ and show diffusive behaviour for filaments, whereas tangent 
vector fluctuations are finite and thus irrelevant for the scaling behaviour of strings. 
In the presence of a potential V = V r + V a + V p , the scaling behaviour of unbound 
segments of a string or filament is governed by the same roughness exponents (provided 
p > 2 for strings and p > 2/3 for filaments [7] ) . 

For unbound strings and filaments, i.e., in the absence of a sufficiently strong 
attractive potential V a , this leads to the scaling form 

Z L = L-*« \z\ e -/ 2 n u (\z\L-S, ML 1 -?) (9) 

in the limit of small \zq\ and |vo|. For strings, the tangent v is an irrelevant scaling 
variable. We introduced exponents Xu characterizing the return probability and 6 U 
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characterizing the segment distribution at z « 0, and a shape function Cl u (y,u) (with 
finite O„(0, 0)) giving the shape of the polymer segment distribution. 

For strings and filaments bound by the attractive potential V a , the longitudinal 
correlation length £|| = 1/\Eq\ gives the characteristic length of unbound segments 
and enters the scaling behaviour, 

Z L = q Xb \z\ e ^ 2 n b Ivllzp-O/C) e i/«n (io) 

with analogous exponents Xb and 9b, which differ from the unbound case in general. 

For a given potential, the two exponents x and 9 are not independent 
as can be seen by using the above scaling forms in the Chapman-Kolmogorov 
relations J J Zl{z\, Vi|z, v)Zl(z, v|zo,Vo) = ^2l(zi,vi|zo,vo) for filaments and 
J z Z l(zi\z) Z l(z|zo) = Z 2 l(z 1 |z ) for strings. This leads to scaling laws 

X = max (dj_/2 + 9/2,0) (strings), X = max (2d± + 39/2, 0) (filaments), (11) 

holding both for Xu, 9 U and Xb, 9b- Exponents x < are not possible because they 
correspond to an unphysical increase of contacts as the length £|| of unbound segments 
increases. If d±/2 + 9/2 < for strings or 2d± +39/2 < for filaments a finite fraction 
of all polymer segments is bound at z — and the main contributions to the z- integrals 
in the Chapman-Kolmogorov relations come from small scales |z| ~ l a leading to x = 
in (JTTJ). 

5. Duality mapping 

Inspecting the Laplace transformed TM equation © and the stationary TM equation 
Q for strings, we observe a formal similarity if we identify s = —E: a short-range 
attractive potential V a (z) cx — <5(z — zq) in the stationary TM equation (j4]) plays the 
role of the initial condition in the Laplace transformed TM equation (j3|) for a potential 

V — V a , i.e., in the absence of the short-range attraction V a - A similar observation 
can be made for the corresponding TM equations and © for filaments where a 
short-range attractive potential V a (z, v) cx — S(z — zo)5(v — vo) in the stationary TM 
equation ([8} plays the role of the initial condition in the Laplace transformed TM 
equation y$ for a potential V ~V a . This is the main idea of the present paper and 
will allow us to establish a duality mapping between the stationary TM equation for 
bound states (characterized by the set of exponents 9b and Xb) in a generic potential 

V = V r + V a + V p and the Laplace transformed TM equation for unbound states 
(characterized by the set of exponents 9 U and Xu) in a potential V — V a — V r + V p 
lacking the short-range attractive part. 

A string in a bound state iPe( z ) fulfils the stationary TM equation (j4j for a 
potential V containing the short-range attraction V a (z) = G8(z — Zq) where we 
consider the limit of small |zo|. We compare the stationary TM equation (Q} with 
the Laplace transformed TM equation (|3]) for Z^~ v "-(z\lq) with s = —E and for 
a potential V — V a without the short-range attraction. If we rewrite 8{z — Zo) = 
5{z — zo)ZY~ Va (z\zo)/ZY~ Va ( z o\ z o) we find that both equations are equivalent and 
solutions have the same normalization if the following two conditions are fulfilled: 

^(z)^M E ZY E Va (z\z Q ) with N E 2 = f[Z v _ E Va {z\z )] 2 (12) 

J Z 

-G- 1 = zY E Va (zo\z Q ) = i>l{z )/M B . (13) 



Letter to the Editor 



5 



These two conditions define the duality mapping for strings between TM equations 
for potentials V and V — V a . 

For filaments we proceed analogously for a bound state V , ]?( z j v ) which fulfils the 
stationary TM equation ([5]) for a potential V containing the short-range attraction 
V a (z, v) = G5(z — zo)<5(v — Vo), where we consider the limit of small |z | and |v |. We 
compare the stationary TM equation (0 with the Laplace transformed TM equation 
(0 for ZY~ Va (z, v|zo, vo) with s = —E and for a potential V — V a without short- 
range attraction. Following analogous steps as outlined for strings above, we find the 
following duality mapping for filaments, 

with M E 2 = f f ^ £ y "(z,v|z ,v )Z i : £ l/ ''(z,-v|z ,vo) (14) 

J z J V 

-G- 1 = z!-^(zo,vo|z ,v ) = ^(z ,v a )/M E , (15) 

relating the TM equations for potentials V and V — V a . This exact mapping can 
be generalized to the more general class of potentials V a = G4>(v)5(z) if we use 
the additional assumption that ZY E Va (zq, vo|zo, 0) ~ 5(y ) is a strongly localized 
function of Vo in the limit Zq « 0. This assumption is justified if the scaling function 
fl a (y,u) is exponentially decaying for u > 1 such that ZY_e a (z , v |z , 0) ss for 
tangents |vo| 3> |zo | 1 ^ 3 ■ Then we can integrate both sides of with a kernel 
f &(vo)ZY~ Va (zo, vo|zo, 0) . . ., which finally leads to a generalized duality mapping 

^(z,v) =A^Zr £ K (z,v|z o ,0) 

-2 / / ZV-V al _ ... n \vV- 



with J\[g — I / ^ B l/ »(z,v|z o ,0)Z! B ^(z,-v|z o ,0) (16) 

J z J V 

-G- 1 = f $(v o )zr s y "(z o ,v o |z o ,0) = [ $(v )Vs(zo,vo)/7V E , (17) 

which is valid in the limit z ~ 0. 

The validity of the duality mappings can be confirmed for a number of potentials 
by direct TM calculations for strings [IT] and filaments [HI [B] . The mappings allow us 
to obtain results for the full potential V by solving the Laplace transformed problem for 
the simpler potential V — V a and give direct information on the partition sums ip^ and 
ZY~ Va and thus the segment distributions. The duality mappings generalize exponent 
relations that have been found previously, as we will show in the following section. 
Furthermore, relations (|T3| . flU} and (fT7|) allow us to determine the transition point, 
i.e., the critical potential strength G c , and the exponent v« describing the divergence 
of the correlation length close to the transition, £|| cx \G — G C | _I/ " . 



6. Exponent relations 

Without working out explicit solutions of the TM equations, we can use the duality 
mapping to derive various exact exponent relations. To derive the exponent relation 
for Xu and Xb for strings we study the limit of small \E\ in ([121) . The scaling form 
(0 for the unbound string determines the s-dependence of the singular part of Z^~ Va 
for small s according to Zj^^ ~ s Xu_1 . For < 1 the singular part is the leading 

order contribution; for Xu > 1 the leading order contribution is finite, ZY~ Va ~ const. 
Using the Chapman-Kolmogorov relation, we find from (I12[) the singular behaviour 
Me ~ |-E| 1_X "/ 2 for x-a < 2 for small \E\ and Me ~ const for x.u > 2. Furthermore, 
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~ (E 1 ^ 6 / 2 for small \E\ according to the scaling form (ITU|) . Equating powers of |J5| 
in (|12l) we arrive at the exponent relation 

J max(2-x«,0) for x u > 1 , . 

Xh = \ X« for X« < 1 (18) 

for strings. For filaments, an analogous analysis of relation f|14jl at small \E\ gives the 
same exponent relation (|18p. For strings, relation ([T5]) agrees with direct calculations 
using the TM equations and also applies in the presence of a long-range power- 
law potential V% ~ w|z| -2 (p = 2), where the exponents \ depend continuously on 
w, as can be checked using the results of Ref. jTD] . Also for filaments, (fT5|) agrees 
with direct TM calculations for potentials V = V a and V = K + V„ [5] IT2]. This 
exponent relation has been formulated in Ref. |13j based on a mapping between the 
renormalization group equations for strings and filaments of different dimensionality. 
An equivalent exponent relation has been confirmed numerically in Ref. [§]. 

In order to derive the corresponding exponent relation for 6 U and 6b for strings 
and filaments, we analyze the scaling behaviour of the Laplace transform ZY~ Va of 
the unbound string or filament for small |z| in (fl~2"f and (|14| . respectively. Using the 
scaling form §§§ for the unbound string or filament we find Z^~ Va ~ |j5|( 1- x«)/C+ 9 i»/ 2 
for Xu > 1 and Z^~ Va ~ |z| 9 "/ 2 for Xu < 1- According to the scaling form (|10[) for 
the bound string or filament we have ipE ~ |z| 8b ^ 2 for small |z|. Equating powers of 
|z| in (TT3|) or (|14[) we arrive at the exponent relation 

9 n + 2(1 - Xu)K for X« > 1 (l9) 
% for X« < 1 1 ' 

which holds for strings with £ = 1/2 and x« = d±/2 + 8 u /2 and filaments with <^ — 3/2 
and Xu = 2dj_ + 36 u /2, according to the scaling laws ([TT|) (note that x„ > for the 
unbound case). For Xb > 0, the same exponent relation can be obtained from a linear 
combination of (p~8j) and the two relations which follow from (fTTj) for the exponent 
pairs Xu, @u and Xb,9b, respectively. Again, it can be checked that relation (|19[) agrees 
with direct TM calculations both for strings [HI \W\ and for filaments [6j [12] [13] . 

Now we address the transition point, transition order, and the correlation length 
exponent by analyzing the dependence of the bound state energy E on the potential 
strength G in relations (fT3"]) and ([Tij]) . Setting E = on the right hand side we find the 
transition point G c . As the singular part of Z^~ Va for small s is ZY~^ ~ s x " _1 , we 
find G c = for Xu < 1; thus, there is no unbinding transition for x„ < 1 and strings 
and filaments are always in a bound state. Expanding around E = for x« > 1 gives 
IG- 1 - G- 1 ! oc {El 1 ^ = ti~ 1/n with 

1/f H = min (xu - 1, 1) for % u > 1 . (20) 

We also used that the linear order dominates the singular contribution to ZY~ Va for 
Xu > 2 such that the transition becomes /srst order with i/m = 1. For 1 < Xu < 2, 
we find i/ii > 1 and a continuous transition. The result (|20[) agrees with those of 
the necklace model [M]. For filaments, relation ([20]) can be generalized for a class 
of tangent-dependent potentials V a = G$a(v)<5(z), satisfying a homogeneity relation 
<1>a(&v) = 6 _a< 1 ) a(v), which has been considered also in Ref. [6j. Performing the 
analogous expansion in (fTTj) we find 



l/z^l = min(x« - 1, 1) for £„ > 1 , where Xu = Xu - d±(l - A)/2 . (21) 
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For this class of potentials there is no transition for x„ < 1, a first order transition for 
Xu > 2 and a continuous transition for 1 < Xu < 2. The result (|2"0)l is recovered for 
A = 1 and $i(v) = <5(v). 

The exponent relations (|18l) and (|2"0)l or (|21l) . together with the scaling law 
allow us to calculate all critical exponents of the unbinding problem if only one 
exponent (xu or 9 U ) of the unbound string or filament in the absence of the short- 
range attractive potential is known. These exponents are often known analytically, or 
can be easily obtained numerically. For V = 0, we have 9 U = #o = for strings and 
filaments. For V = V r and d± = 1, we can make use of another exponent relation, 
Xu = Xr = 1 + C [E]; which is also valid for both strings and filaments. 

7. Conclusions 

In conclusion we derived a duality mapping between bound and unbound states of 
one-dimensional strings and filaments. This mapping allows us to determine the 
transition point and the order of unbinding and desorption transitions of strings and 
filaments. We derived exponent relations for the return probability exponents x> 
the segment distribution exponents 9 and the correlation length exponent v\\ from 
the mapping. These relations allow us to determine all critical exponents related to 
the unbinding and desorption transitions of both filaments and strings from a single 
exponent characterizing the unbound string or filament. 
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